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ABSTRACT 

We use an analytic model to investigate the theoretical uncertainty on the thermal Sunyaev- 
Zel'dovich (SZ) power spectrum due to astrophysical uncertainties in the thermal structure of the 
intracluster medium. Our model accounts for star formation and energy feedback (from supernovae 
and active galactic nuclei) as well as radially dependent non-thermal pressure support due to random 
gas motions, the latter calibrated by recent hydrodynamical simulations. We compare the model 
against X-ray observations of low redshift clusters, finding excellent agreement with observed pressure 
profiles. Varying the levels of feedback and non-thermal pressure support can significantly change 
both the amplitude and shape of the thermal SZ power spectrum. Increasing the feedback suppresses 
power at small angular scales, shifting the peak of the power spectrum to lower i. On the other hand, 
increasing the non-thermal pressure support has the opposite effect, significantly reducing power at 
large angular scales. In general, including non-thermal pressure at the level measured in simulations 
has a large effect on the power spectrum, reducing the amplitude by 50% at angular scales of a few 
arcminutes compared to a model without a non-thermal component. Our results demonstrate that 
measurements of the shape of the power spectrum can reveal useful information on important physical 
processes in groups and clusters, especially at high-redshift where there exists little observational data. 
Comparing with the recent South Pole Telescope measurements of the small-scale cosmic microwave 
background power spectrum, we find our model reduces the tension between the values of erg measured 
from the SZ power spectrum and from cluster abundances. 

Subject headings: cosmology: dark matter — galaxies: clusters: general — intergalactic medium 
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1. INTRODUCTION 

The Sunyaev-Zel'dovich (SZ) effect has long been rec- 
ognized as a powerful tool for probing the physics of the 
intracluster medium (ICM), large-scale structure forma- 
tion and the dark energy equation of state (Birkinshaw 
1999; Carlstrom et al. 2002). Experiments such as the 
South Pole Telescope (Ruhl et al. 2004), the Atacama 
Cosmology Telescope (Kosowsky 2003) and Planck^ are 
currently surveying the microwave sky with the goal of 
identifying clusters via their SZ signature on the cosmic 
microwave background (CMB), thus constructing large 
catalogs of galaxy clusters that are uniformly selected by 
SZ flux (Vanderlinde et al. 2010). 

The SZ effect can also be detected as a secondary 
anisotropy in the CMB temperature power spectrum, 
appearing as "excess power" (over the predicted pri- 
mary anisotropy signal) on angular scales of a few arc- 
minutes. _ The ensemble averaged power spectrum am- 
plitude Ci has an extremely sensitive dependence on trg, 
oc (T|(f2f,/i)^ (Komatsu & Seljak 2002), where erg is the 
power spectrum normalization on scales of 8/i~^ Mpc. 
The SZ angular power spectrum thus represents a ro- 
bust observable with which competitive constraints on 
CTg can be obtained. This would in turn enable tighter 
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constraints to be placed on the dark energy equation of 
state parameter w by breaking the degeneracy with erg in 
constraints derived from the primary CMB signal alone. 

A number of experiments have reported detections of 
power in excess of the primary CMB power spectrum 
at small angular scales and thus been able to place up- 
per limits on the amplitude of the SZ signal (Dawson 
et al. 2006; Friedman & QUaD Collaboration 2009; Re- 
ichardt et al. 2009a, b; Sievers et al. 2009; Sharp et al. 
2010; Fowler et al. 2010). The South Pole Telescope 
(SPT, Lueker et al. 2010) has recently reported mea- 
surements of the CMB power spectrum at 150 GHz for 
angular scales in the range 2000 < £ < 10000. By 
combining SPT survey maps at 150 and 220 GHz to 
minimize astrophysical foreground signals, Lueker et al. 
(2010) were able to isolate and detect SZ power (ki- 
netic plus thermal) at 2.6a. The measured amplitude 
at ^ = 3000 was 4.2 ± l.hfxK^, significantly below that 
predicted by halo model calculations (Komatsu & Sel- 
jak 2002) or simulations (White et al. 2002; Shaw et al. 
2009; Sehgal et al. 2010), assuming WMAP7 cosmolog- 
ical parameters. The significantly lower-than-predicted 
signal could be explained by a lower value of erg. Us- 
ing the SZ power spectrum predicted by simulations and 
assuming no modeling uncertainty, Lueker et al. (2010) 
measured ug = 0.746 ± 0.017. However, this result is 
in tension with other probes of CTg; for example, re- 
cent X-ray and optical measurements of cluster abun- 
dances suggest CTg = 0.82 ± 0.05 (Mantz et al. 2010), 
CTg(r2ji//0.25)° '^^ = 0.813 ±0.013 (Vikhlinin et al. 2009), 
and CTg = 0.807 ± 0.02 (Rozo et al. 2010). 

This disagreement can be resolved if current analytic 
models and simulations over-predict the level of SZ power 
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by a factor of « 2. SZ power spectrum calculations have 
two sources of uncertainties: the amplitude of the halo 
mass function, and in modeling the radial electron pres- 
sure profile of the ICM. Recent studies of the mass and 
redshift distribution of halos in N-body simulations im- 
ply that the mass function is known to 5%-10% accuracy 
for the currently allowed wCDM cosmology (Tinker et al. 
2008; Bhattacharya et al. 2010a) and so is not the dom- 
inant theoretical uncertainty (although, as pointed out 
by Stanck et al. 2009, the impact of baryonic physics on 
halo masses is still somewhat unclear). 

The main difficulty in calculating the thermal SZ (tSZ) 
power spectrum is predicting the thermal pressure pro- 
files of groups and clusters over a wide range of mass and 
redshift. Komatsu & Scljak (2002) and Holder (2002) 
show that low-mass (Af < 2 x 10" h^^Mo) and high- 
rcdshift (z > 0.6) objects both make a significant con- 
tribution to the signal at angular scales of a few arcmin- 
utes {t f« 4000), the scales at which current small-scale 
anisotropy experiments such as SPT are most sensitive. 
While Chandra and XMM-Newton have enabled high- 
mass, low- redshift clusters to be studied extensively over 
the last decade, lower mass and, in particular, high- 
redshift objects have not been extensively studied. Fur- 
thermore, the SZ power spectrum is sensitive to the ther- 
mal pressure of the ICM out to cluster radii several times 
larger than is typically probed by X-ray observations. 
Hence, there are few observational constraints that can 
be used to aid predictions of the tSZ power spectrum. 

From a theoretical perspective, full cosmological hy- 
drodynamical simulations have only recently begun to 
systematically explore the effects of sub-grid baryonic 
processes, such as radiative cooling, star formation, feed- 
back mechanisms, cosmic rays, thermal conduction and 
magnetic fields, on the thermal structure of the ICM (Na- 
gai et al. 2007a; Pfrommer et al. 2007; Dolag & Stasyszyn 
2009; Sijacki et al. 2007, 2008; Battagha et al. 2010). Fur- 
thermore, simulations require high spatial resolution in 
order to effectively implement some of these processes. 
However, large box sizes are also required in order to 
adequately sample the halo mass function at group and 
cluster scales to enable measurements of the tSZ power 
spectrum. Currently, the computational expense of run- 
ning large box, high-resolution hydrodynamical simula- 
tions is prohibitive to investigating the level of theoret- 
ical uncertainty on the power spectrum as well as the 
dependence on cosmological parameters. 

The principal aim of this work is to use analytic mod- 
els to investigate variations in the predicted tSZ power 
spectrum caused by uncertainties in the thermal struc- 
ture of the ICM. Specifically, we study the impact of 
energy feedback, non-thermal pressure support and halo 
concentration on the shape and amplitude of the power 
spectrum. This is achieved by combining an analytic 
model for the ICM with the halo mass function to rapidly 
calculate the power spectrum for different model param- 
eters. Our model assumes that gas resides in hydrostatic 
equilibrium with a polytropic equation of state, and ac- 
counts for star-formation as well as feedback from super- 
novae and active galactic nuclei (AGNs) and non-thermal 
pressure support driven by random gas motions and tur- 
bulence in the ICM. We calibrate our model parame- 
ters by comparing against X-ray observations of mas- 



sive, low-z clusters. We find that including non-thermal 
pressure support at the level measured in state-of-the- 
art hydrodynamical simulations significantly reduces the 
amplitude of the predicted tSZ power spectrum, thus re- 
ducing the tension between the as inferred from the SPT 
observations and cluster abundance measurements. 

The paper is organized as follows. In Section 2, we 
describe our model of the ICM. In Section 3, we explore 
the model parameter space by comparing model scaling 
relations and radial profiles against recent low-redshift 
X-ray observations of groups and clusters. In Section 4, 
we explore the uncertainty on the tSZ power spectriim 
sourced by the underlying range in ICM model param- 
eters. We compare our fiducial power spectrum model 
with other recent simulations, and discuss our results in 
the context of the recent SPT observations. 

Throughout this paper we assume a fiducial, spa- 
tially flat, ACDM cosmological model consistent with 
the WMAP7 best-fit cosmological parameters, namely 
Ho = 71 km s-^Mpc-i, = 0.264, = 0.044, 
Oa = 0.736, Hs = 0.96 and erg = 0.8. 

2. THEORETICAL MODELS 

2.1. Thermal SZ Power Spectrum 

The tSZ effect is a distortion of the CMB caused by 
inverse Compton scattering of CMB photons off electrons 
in the high temperature plasma within galaxy clusters. 
To first order, the temperature change at frequency v 
of the CMB is given by AT/Tcmb{xu) = f{xv)y, where 
J{xv) = a;i.(coth(a;,./2) - 4), x„ = hv/kBTcuB, and y is 
the dimensionless Compton-y parameter 

y = (^) I nS)Te{l)dl , (1) 

where the integral is along the line of sight, Tcmb the 
CMB temperature, and rig and Tg are the number density 
and electron temperature of the ICM, respectively. 

The tSZ power spectrum can be calculated by simply 
summing up the squared Fourier-space SZ profiles of all 
clusters: 

C, = f{x.fj dz^ I d]nM^^^^y\M,z,i) (2) 

where V{z) is the comoving volume per steradian and 
n{M, z) is the number density of objects of mass M at 
redshift z. For the latter we use the fitting function of 
Tinker et al. (2008). y{M, z, £) is the Fourier transform of 
the projected SZ profile for a cluster of mass M and red- 
shift z. This can be calculated assuming spherical sym- 
metry using (Bracewell 2000; Komatsu & Seljak 2002): 

.^ 47rrc /"^ , Or./,. ,sin(fo/4) 
y{M,z,£) = -^ dxx^Py{M,z,x)—y^ , 

' (3) 

where x = rjrc, ic = DA{,z)/rc., rc is a characteristic 
scale radius of the profile, Da{z) is the angular diam- 
eter distance to redshift z and py{M,z,x) is the three- 
dimensional SZ profile, which is related simply to the gas 
thermal electron pressure profile, 

Py{M,z,x) = ^Pe{M,z,x) . (4) 
Note that while this calculation only accounts for the 
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one-halo contribution to CV, Komatsu & Kitayama 
(1999) demonstrated that, at the angular scales being 
probed by the current generation of small-scale CMB 
experiments (i.e., t > 1000), the two-halo (or clustered) 
contribution to the tSZ power spectrum is nearly 2 or- 
ders of magnitude smaller than the one-halo term and 
decreases with increasing £. We thus neglect the two- 
halo contribution in our calculation. 

Groups and clusters over a wide mass (10^^ < 
M/M© < 10^''^) and redshift (0 < z < 3) range con- 
tribute to the tSZ power spectrum. For example, more 
than half the power at £ = 3000 comes from clusters 
at redshift greater than one (Komatsu & Seljak 2002). 
Likewise, at the same angular scale, half of the power 
is predicted to be sourced by objects of mass less than 
2 X IO^^Mq. While the mass and redshift contribution 
depends on the details of the underlying gas physics in- 
corporated in models and simulations, it is clear that 
much of the signal comes from objects for which there 
exists little direct observational data. 

2.2. Cluster Model 

Our model for the density and temperature structure of 
the ICM is based upon the model of Ostrikcr ct al. (2005, 
henceforth O05) and Bode et al. (2009, henceforth B09), 
with an important modification that allows for a radially- 
dependent non-thermal pressure component (e.g., turbu- 
lence or bulk flows). In this section, we briefly review the 
model and describe how non-thermal pressure support is 
implemented; we refer the reader to O05 and B09 for a 
more detailed description of the original model. 

2.2.1. Dark Matter Halo Structure 

The model assumes that the ICM initially follows 
the density and temperature of the host dark matter 
halo, but rapidly rearranges into hydrostatic equilibrium 
within the potential well of the host dark matter halo 
with a polytropic equation of state. 

Dark matter halo properties are determined by a 
Navarro-Frenk- White (NFW) density profile (Navarro 
et al. 1996, 1997) of the form. 



sity; 



PDM{r) 



(5) 



where x = r/rs, Vg is the NFW scale radius, and ps 
is a normalization constant. Halo concentration is de- 
fined as c = Rvir/rs- Numerous N-body simulations 
have demonstrated there to be a power-law scaling be- 
tween halo mass, redshift and concentration. We adopt 
the halo mass-concentration relation measured by DuSy 
et al. (2008) from halos identified in a large N-body sim- 
ulation over the redshift range < ^ < 2, 



c(M, z) = 7.S5Ac 



2 X 1012 /i-iMm 



-0.081 



-0.71 



(6) 

where Ac is an arbitrary normalization factor that we 
will later vary to investigate the impact of the mass- 
concentration normalization on the tSZ power spectrum 
through our model. Ac = 1 gives the mass-concentration 
relation of Duffy et al. (2008). 
Halo masses are defined in terms of the virial overden- 
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where Ac = IStt"^ + 82{nM{z) - 1) - 39(17m(-2) - 1)^ is 
the virial overdensity given in Bryan & Norman (1998) 

and Pc{z) is the critical density at redshift z. When com- 
paring with observations, we also use the mass definition 
-^-^500 = (4/3)7ri?5QQ.500/;(.: where i?5oo is the radius out 
to which current X-ray observations can reliably measure 
the gas density and temperature profiles. 

2.2.2. Star Formation 

The gas is initially treated as a tracer of negligible 
mass of the dark matter gravitational potential, i.e., it 
does not contribute itself to the gravitational potential of 
the cluster and the gas density and temperature follows 
that of the dark matter. 

We assume that some fraction of the gas has radia- 
tively cooled and formed stars. To determine the stellar 
mass of a cluster we use the observed scaling relation 
of Giodini et al. (2009), who measured the ratio of stel- 
lar mass to total mass for 91 X-ray selected groups and 
low-mass clusters in the COSMOS survey (Scoville et al. 
2007). This was supplemented by the 27 low- redshift. 
X-ray selected clusters analyzed by (Lin et al. 2003), re- 
sulting in an overall sample that enconipassed a wide 
range of both mass (10^3 < M500 < lO^^) and redshift 
(0 < z < 1). They found that the stellar mass fraction 
within i?5oo, /*(< i?50o) = M*(< i?50o)/-^^500, followed 
a mean empirical relation; 



M. 



500 



3 X 1014 M, 



-0.37±0.04 



/* = (2.58 ±0.05) X 10" 

We assume that the stellar mass fraction given above 
also holds within _Rvir, i-o., /* = M»(< i?5oo)/M5oo = 
M»(< i?vir)/Mvir. Following O05 and B09, we deter- 
mine the redshift evolution of the stellar mass fraction 
adopting the "fossil" model of Nagamine et al. (2006). In 
this model the star-formation rate is given by a delayed 
exponential, with a decay time of 1.5 Gyr for bulge pop- 
ulations and 4.5 Gyr for disk populations. This model 
predicts only a small (« 15%) decrease between z = 
and 1 in the normalization of this relation, which is con- 
sistent with the results of Giodini et al. (2009). 

The initial total gas mass within the virial radius is 
thus 



(/fc-./*)Mv 



and the initial total energy of the gas is 



Eg,i = fb 



2tt 



PDM{r)3al,M{r)r'^dr + 



r 



(9) 



<^ir)——dr 
dr 



(10) 

where fb = flb/^M, $(r) is the gravitational potential 
and (JDMir) is the one- dimensional dark matter veloc- 
ity dispersion (see, for example. Equation (13) of Lokas 
& Mamon 2001, and we assume isotropic dark mat- 
ter orbits). Thus, the initial total energy of the gas 
is simply the sum of the kinetic and potential energy 
of the dark matter halo scaled by the cosmic baryon 
fraction. We assume that the gas within a radius 
has cooled and formed stars, where r* is determined by 
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fbMDMi< '^^,) = /*Mvir and so omit the contribution 
to the total gas energy within this radius {r^, ^ 0.24 
Mpc for a 2 X 10^'* H'^Mq cluster at z = 0.5, decreasing 
weakly toward higher redshift /lower mass) . We note that 
the final gas mass within i?vir predicted by our model is 
not necessarily equal to Mg^i, as the gas may expand or 
contract when hydrostatic equilibrium is enforced (see 
Section 2.2.4). 

2.2.3. Hydrostatic Equilibrium 

The ICM should rapidly rearrange itself into hydro- 
static equilibrium, satisfying 



dPtotjr) 
dr 



d^{r) 
dr 



(11) 



where Pg{r) is the gas density at radius r from the clus- 
ter center and Ptot is the total gas pressure, Ptotir) = 
Pth{r) + Pntir)- The total gas pressure is therefore a 
combination of the thermal and non-thermal pressure 
components, with the latter primarily due to random 
gas motions and turbulence in the ICM. We discuss how 
we implement a non-thermal pressure component in the 
following section. 

Hydrodynamical simulations indicate that the total 
pressure adheres more closely to a polytropic equation 
of state than the thermal pressure, i.e., Ptot oc p^. Fig- 
ure 1 shows the pressure-density phase diagram obtained 
for 16 clusters simulated using the Eulerian hydrodynam- 
ics ART code (Kravtsov et al. 2002), including radiative 
cooling and star formation (Nagai et al. 2007a). The to- 
tal pressure (blue line) is the sum of the thermal (red line) 
and non-thermal components in radial bins around the 
cluster potential minimum, averaged over the entire sam- 
ple. The shaded regions denote the standard deviation 
within the sample around the mean. The non-thermal 
pressure is measured from the radial velocity dispersion 
of the gas within each shell (Lau et al. 2009). The black 
dashed line represents F = 1.2. 

It is clear from this figure that using the total pressure 
results in a constant adiabatic index of 1.2 over more 
than four decades of gas density, whereas the thermal 
pressure implies a density (and hence radially) depen- 
dent F. The same results are obtained for non-radiative 
simulations without gas cooling and star-formation. We 
therefore assume that the total gas pressure has a poly- 
tropic equation of state, Ptot = Po{p / Po)^ , where we set 
F = 1 + 1/rt = 1.2 (n = 5, where n is the polytropic in- 
dex) and po and Pq are the central density and pressure 
of the gas. In Section 2.2.5 we describe our method for 
including non-thermal pressure in the model and demon- 
strate that this causes the ratio Pth/p^ to vary with clus- 
ter radius. 

The final total pressure and density of the gas in our 
model are given by 



Ptot{r)^PoO{rr+^ 
Pgir)^po9{ry\ 

where 9{r) is the polytropic variable 

t M) 

and <l>o is the central potential of the cluster. 



(12) 
(13) 



(14) 
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Fig. 1. — Total (blue line) and thermal (red) pressure as a func- 
tion of gas number density, averaged over 16 clusters simulated by 
Nagai et al. (2007a) at z = 0. The shaded regions show the stan- 
dard deviation of the cluster sample around the mean. The solid 
black line shoves the best fit line to the total pressure and density 
with a slope (adiabatic index) of F = 1.2. The lower panel shows 
the fractional deviation of the pressure around this line. 

2.2.4. Energy Feedback 

Once Pq and po are determined, the density and pres- 
sure profile of the gas is fully specified. These variables 
are determined by applying two constraints. First, the 
total final energy of the gas that is initially within i?vir 
must obey 



DM] 



(15) 



This condition states that the final energy Egj be equal 
to the initial energy, Eg i (Equation 10), plus energy 
added to the gas via dynamical process or from en- 
ergy feedback (discussed below), plus the work done by 
the gas as it expands or contracts relative to its ini- 
tial state. We define Rf as the radius within which 
the final gas mass, following polytropic rearrangement, 
is equal to Mg^i (Equation 9). The work done by the gas 
is then given hy AEp = (47r/3)(i?3.^ - Rj)P,. Hence, 
if i?/ > i?vir the gas has expanded, doing work in 
the process. Pg is the surface pressure at i?vir, given 

by Ps = fbPDM{Rvir) = fbO-DMiRvi':)PDMiRvh)- The 

surface pressure of the gas at the virial radius is thus 
equal to that of the dark matter multiplied by the cos- 
mic baryon fraction. The second constraint used to 
solve the model is that the total pressure of the gas at 
Rf he equal to the gas pressure at the virial radius of 
the cluster before polytropic rearrangement, i.e., that 

Ptot{Rf) = hPoAliRvir)- 

One fundamental component of the model is the in- 
clusion of heating of the ICM via non-gravitational pro- 
cesses. This can occur via two mechanisms; (1) through 
energy transfer from the dark matter to gas during ma- 
jor mergers, and (2) via energy feedback from supernovae 
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and AGN outflows. 

To account for the former, Bode et al. (2009) intro- 
duced the parameter cdm which controls the fraction of 
the total dark matter energy, \Edm\^ transferred to the 
ICM during mergers (Pcarce et al. 1994; Rasia ct al. 2004; 
Lin ct al. 2006; McCarthy et al. 2007). The total dark 
matter energy is simply the sum of the total dark mat- 
ter kinetic and potential energy (calculated similarly to 
Equation 10, but setting r* = and fh = 1). Prom the 
hydrodynamical simulations of McCarthy et al. (2007), 
B09 suggest EDM = 0.05. 

Energy feedback from supcrnovae or AGNs is assumed, 
as a first approximation, to be proportional to the to- 
tal stellar mass of a cluster and is determined by the 
parameter cf. The approximate value of this parame- 
ter is difficult to determine from hydrodynamical simu- 
lations as methods of implementing AGN feedback have 
only recently started to be investigated in earnest (Si- 
jacki et al. 2007; Di Matteo et al. 2008; Booth & Schaye 
2009; Teyssier et al. 2010). B09 demonstrated that the 
feedback parameter can be calibrated by comparing the 
model to low-redshift X-ray scaling relations, finding 
values in the range < Cf < 1.2 X 10"^ depending 
on assumptions of the stellar mass of clusters and its 
mass dependence. Recently, Battaglia et al. (2010) in- 
vestigated the impact of AGN heating in hydrodynam- 
ical simulations, finding an effective feedback efficiency 
(which is roughly equivalent to our ef ) of 5 x 10~^ for an 
-^500 = 6.8 X 10^^ h~^MQ cluster over the duration of 
the simulation. 

We note that, in terms of predicting the tSZ power 
spectrum, the effect of cf and cdm is somewhat degen- 
erate. Following B09, we take cdm = 0.05 as our fidu- 
cial parameter. In the following section we compare our 
model using different values of the feedback parameters 
et and eDM with direct X-ray observations of clusters. 
In Section 4, we investigate the impact of varying both 
parameters on the SZ power spectrum. 



2.2.5. Non-thermal Pressure 

Hydrodynamical simulations have demonstrated that 
a significant fraction of the total energy of the ICM is 
contained within random gas motions, which provide a 
significant non-thermal contribution to the total pres- 
sure support (Evrard 1990; Rasia et al. 2004; Kay et al. 
2004; Dolag et al. 2005; Lau et al. 2009). This 'kinetic 
support is largely sourced by infalling and merging struc- 
tures, which can further generate turbulent gas motions 
at the boundary between the bulk flow and the ther- 
malized ICM (Dolag et al. 2005; Kim 2007; Vazza et al. 
2009). While the level of non-thermal pressure is typi- 
cally found to be small in the central regions of clusters, 
it increases steadily with radius, becoming a significant 
fraction of the total pressure at R500 (Lau et al. 2009). 
The simulation comparison study of Prenk et al. (1999) 
demonstrated that different hydrodynamics codes pre- 
dict a similar ratio of kinetic to thermal gas energy in 
clusters. 

To account for non-thermal pressure support in our 
model, we split the total pressure into non-thermal and 

thermal components, Ptot{f) = Pth{'>') + P-ntif) and set 
the non-thermal pressure fraction to be a power law with 




r/Rsoo 

Fig. 2. — (Upper) The ratio of non-thermal to total pressure as a 
function of radius for the 16 simulated clusters of Lau ct al. (2009) 
(gray lines). The thick solid line is the best fit from Equation 16 
and represents the fiducial non-thermal pressure fraction assumed 
in our model. (Lower) The thermal polytropic index V^f^ (sec 
text) predicted by our gas model as a function of radius. The 
dotted, solid and dot-dashed black line represents ao = 0, 0.18 and 
0.3, respectively. The dashed line represents the fiducial model 
ao = 0.18, but with no star-formation or feedback. 

cluster-centric radius; 

^(z) = a(z) , (16) 

where n„t determines the radial dependence of the non- 
thermal pressure fraction. 

We assume that the non-thermal pressure fraction 
varies with redshift according to a{z) = aof(z), where 
Q:o is the mean ratio of non-thermal to total presstire at 
i?50o at z = 0, and f{z) is a monotonically increasing 
or decreasing function of redshift. Since Pnt/Ptot < 1, 
Equation 16 is valid for a{z) < (iJmax/.Rsoo)"""*- As 
we describe below, the outermost radius of the pressure 
profiles in otir model is set to 4i?5oo, which limits the 
maximum value to a{z) < 4~""*. We adopt the follow- 
ing form for the redshift evolution, 

f{z) = min[(l + zf, ifmax - 1) tanh(/3^) + 1] , (17) 

where fmax = 4~"''*/q!o and /3 is a free parameter that 
determines the evolution rate. This form ensures that, at 
low redshift {z <1), a{z) has a simple power-law depen- 
dence on (1 -|- 2;), whereas at high redshift (or for large 
values of ao) it smoothly asymptotes toward the maxi- 
mum value of 4^""*. We note that l3 can be positive or 
negative. For the former, the value of f3 controls the red- 
shift of transition from the power law to the asymptotic 
behavior of a{z). For negative values of /?, a{z) declines 
toward zero with increasing redshift. 

We use the hydrodynamical simulations of Nagai et al. 
(2007a) to find the best-fit values of ao and Unt- The 
upper panel of Figure 2 shows the radial dependence of 
the ratio of the non-thermal to total gas pressure frac- 
tion, Pnt/Ptot, of 16 simulated clusters a.t z = (the 
same set of simulations shown in Figure 1). Following 
Lau et al. (2009), the non-thermal pressure component 
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(due to random gas motions) was measured by determin- 
ing the velocity dispersion of gas cells in radial shells in 
the rest frame of the cluster. It is clear that the mean 
non-thermal pressure fraction increases rapidly with ra- 
dius, from < 10% at r = 0.1i?5oo to greater than 30% 
at r = 2i?5oo. The thick-solid line indicates the best-fit 
model with ao = 0.18 ± 0.06 and n„t = 0.8 ± 0.25 (re- 
sulting in a maximum value of a{z) = 0.33). We find 
no evidence of mass dependence on any of the param- 
eters. Furthermore, Lau ct al. (2009) obtained similar 
results for their non-radiative simulations (i.e., includ- 
ing no cooling or star formation) oiitside the central re- 
gions of the clusters. This indicates that the level of non- 
thermal pressure is not sensitive to these processes. We 
note, however, that the preheating simulations of Stanek 
et al. (2010) do suggest a lower level of kinetic energy in 
bulk motions than in simulations without preheating. 

We determine the redshift dependence by examining 
a high- redshift output of the simulation. At z = 1, we 
find that the best-fit ao increases to 0.26, implying that 
^ « 0.5. A greater fraction of non-thermal pressure to- 
wards high redshift is expected due to the increased rate 
of merger activity in the ACDM model. Stanek et al. 
(2010) also observed an increasing level of kinetic support 
with redshift in their 'gravity-only' smoothed-particle hy- 
drodynamics simulation. However, the results of their 
preheating simulation imply a decreasing (/3 < 0) level 
of non-thermal pressure toward high redshift. The spe- 
cific values of ao and ,5 prcidicted by simulations thiis 
appear to be somewhat dependent on the baryonic pro- 
cesses included. We therefore vary these parameters to 
investigate their impact on the tSZ power spectrum. Fur- 
ther study is required to determine the dependence of 
non-thermal pressure to cooling and heating effects in 
hydrodynamical simulations. 

We find that n,,/ decreases by approximately 15% at 
z = 1 compared to the best-fit z = value; however, as 
this change is less than the 0.25 standard deviation on 
the redshift zero value, we henceforth hold n„t constant 
with redshift. Our fiducial model is therefore ao = 0.18, 
nnt = 0.8 and /3 = 0.5. We investigate the impact of 
varying the overall normalization of the profile ao and /? 
on the electron pressure profiles of clusters and the tSZ 
power spectrum, while fixing n„( to the fiducial value for 
the remainder of the paper. In a future work, we intend 
to investigate in more detail the mass and redshift de- 
pendence of non-thermal pressure profiles in much larger 
samples of simulated groups and clusters. 

The radial thermal pressure profile given by our model 

is 



Pth{r) = Podir) 



n+l 



1 - a{z) 



R. 



■500 



0.8 



(18) 



The gas density profile remains as in Equation 13. The 
radial temperature profile is then given by 



kT{r) = iJLmp—e{r) 
Po 



1 - a{z) 



500 



0.8 



(19) 



where /x is the mean molecular weight and rup is the 
proton mass. 

One consequence of including of non-thermal pres- 
sure support in our model is that the adiabatic index 



of the thermal component of the gas, given by Tth = 
log(Pt/i/Po)/log(/Og/po), varies as a function of cluster 
radius. X-ray observations of clusters (Vikhlinin et al. 
2006; De Grandi & Molendi 2002) as well as hydrody- 
namical simulations (Battaglia et al. 2010) indicate that 
the ICM thermal pressure and density profiles do not ad- 
here to a polytrope with constant index. The lower panel 
in Figure 2 shows Tth as a function of radius produced 
by our model for a M500 = 3 x 10^* h~'^MQ cluster. 
The solid lines shows the results for ao = 0.18, the dot- 
dashed line for ao = 0.30 and the dotted line for ao = 0. 
The dashed line gives the results for ao = 0.18, but as- 
suming Cf = cdm = Af* = 0. In the zero non-thermal 
pressure case, we obtain Tth = F = 1.2, as expected. 
Setting a = 0.18 (0.3) results in a w 5 (8)% increase in 
Tth, which also increases with radius. Turning off star 
formation and feedback in our model reduces this radial 
dependence slightly. 

3. COMPARISON WITH X-RAY OBSERVATIONS 

A fundamental test of our model is to compare with 

the global properties and radial profiles of clusters ob- 
tained from X-ray observations. The main aim is to 
evaluate the range in which the two principal parame- 
ters in this model, the feedback parameter cf and the 
non- thermal pressure support parameter ao, reproduce 
observed cluster properties. We focus specifically on the 
M^OQ — fg relation measured from the cluster samples of 
Vikhlinin ct al. (2006) and Sun et al. (2009) (where the 
gas fraction fg is defined as the gas mass divided by total 
mass within iisoo), and the 'miiversar electron pressure 
profile measured by Arnaud et al. (2010) from the REX- 
CESS cluster sample (Bohringer et al. 2007; Pratt et al. 
2009). The latter provides an important test as the ra- 
dial thermal pressure profile is the cluster property that 
determines the SZ power spectrum. We have chosen to 
use the M500 — fg relation as we find the gas fraction 
in our model to be particularly sensitive to the precise 
values of Cf. 

3.1. Hydrostatic Mass Estimates 

One of the principal techniques for measuring cluster 
masses is to measure the radial gas density and tem- 
perature profiles and solve the equation of hydrostatic 
equilibrium to derive the total mass profile, 



M{< r) 



Gpg dr 



(20) 



(Sarazin 1986; Evrard et al. 1996). However, direct ob- 
servations of clusters only currently probe the thermal 
pressure profile and so miss the significant contribution 
of non-thermal pressure in the ICM. This can result in a 
systematic underestimation of M500 by 10% — 20% (e.g., 
Rasia et al. 2004; Nagai et al. 2007b; Piffaretti & Val- 
darnini 2008; Meneghetti et al. 2010). Hence one must 
take care when comparing simulation or model scaling 
relations with observations that utilize this technique. 

Vikhlinin et al. (2006) and Sun et al. (2009) use this 
method to measure the masses of the groups and clusters 
in their sample. Arnaud et al. (2010) determine the mass 
dependence of the amplitude of the pressure profiles of 
the clusters in their sample using the Yx — Afsoo relation 
calibrated by Vikhlinin et al. (2006). Therefore, in order 
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Fig. 3. — Ratio of hydrostatic to true mass measured within 
spheres of radius r/Rsoo for a cluster of M500 = 3 X 10"'^'' h~^MQ. 
We vary the values of the non-thermal pressure support parameter, 
< «o < 0.3 (top to bottom) in steps of 0.06. The solid line 
represents the fiducial value ao = 0.18. The blue point represents 
the mean M|^q^/M5oo measured from the simulated clusters of 
Lau et al. (2009), with the error bar representing the error on the 
mean. 

to compare our model with the observational data, we 
use the hydrostatic mass predicted by our model 
rather than the true M^qq. 

Figure 3 shows the ratio of hydrostatic to true mass 
obtained from our model, varying the value of the non- 
thermal pressure parameter, ag from zero (top line) to 
0.3 (bottom line) in steps of 0.06. The thick solid line 
represents our fiducial value of a = 0.18. The blue point 
represents the mean / Mc^qq measured from the Lau 

et al. (2009) simulated clusters, with the error bar rep- 
resenting the error on the mean. 

For our fiducial model, underestimates the true 

mass M500 by a factor of 0.87 at r = R500, increasing 
to 0.75 at r = 2i?5oo- Setting ao — 0.3 increases this 
to 0.78 and 0.59 at these radii, respectively. We find 
that changing the feedback parameter ef has only a small 
effect on the hydrostatic mass estimate. 

Note that the dependence of M||,q^ on ao presents a 
complication in comparing the model with observations 
as holding M^q^ constant requires varying the true mass 
-^^500 as we vary the model parameters. This makes it 
difficult to determine whether changes in predicted prop- 
erties are due to changes in the true mass or model pa- 
rameters. To overcome this, in the remainder of this sec- 
tion we fix M|J|q^ = O.87M500 (the ratio obtained from 
our fiducial value of ao) for all model realizations. 

3.2. M|J^^ - fg relation 

In Figure 4, we compare the — fg relation pre- 

dicted by our model for different values of ef against the 
X-ray observations. The red points are from Vikhlinin 
et al. (2006), whereas the blue points are from Sun 
et al. (2009). The solid lines represent the results from 
our ICM gas model with (from top to bottom) ef = 
10~^, 10~^, 5 X 10~^, and 10~^. For comparison, we also 
plot results for ef = eDM = (dashed line) and a maximal 



feedback model with cdm = 0.1, et = 10~^ (dot-dashed 
line). In this plot, ao is kept fixed at the fiducial value 
of 0.18 (note that as we assume a polytropic equation 
of state between gas density and total pressure, chang- 
ing the non-thermal pressure fraction does not change 
fg). The dotted line indicates the cosmic baryon frac- 
tion, ilb/^M- We have verified that our model predicts 
fg w flb/flM in the absence of star formation and feed- 
back. 

It is clear that for high-mass clusters varying the feed- 
back parameter ef (solid lines) produces little change in 
the gas fraction. The energy added to the ICM is small 
compared to the total binding energy for these clusters 
and thus does not disrupt the gas distribution signifi- 
cantly. However, the impact of energy feedback becomes 
increasingly important towards lower mass. As cluster 
mass decreases the feedback energy becomes a larger 
fraction of the total binding energy. In these systems, 
feedback has the effect of 'inflating' the gas distribution, 
thus reducing the baryon fraction within a fixed radius. 

Setting eDM = (dashed line) demonstrates that the 
dark matter energy transfer parameter has a constant 
eff'ect with mass, reducing fg by 0.01 — 0.02 across the 
entire mass range. This is because the energy transferred 
to the gas from the dark matter is proportional to the 
binding energy of the cluster, and so has a more signifi- 
cant effect on the high mass clusters than ef. 

There is generally a large amount of intrinsic scatter 
in the observed sample, especially at the high-mass end 
where the error bars on individual estimates of fg are 
smaller. Although the Vikhlinin et al. (2006) clusters 
were selected for their relaxed morphology, the large scat- 
ter may indicate a range of merger histories within the 
sample. It may also reflect the large (« 50%) intrinsic 
scatter observed in the mean stellar mass - total mass 
relation (Giodini et al. 2009; Gonzalez et al. 2007). Nev- 
ertheless, the upper and lower model proflles plotted in 
Figure 4 bracket the observed data points. We therefore 
assume this range in ef (0-10~^) and eDM (0-0.1) when 
evaluating the modeling uncertainty on the tSZ power 
spectrum. 

3.3. Electron Pressure Profiles 

The two key components required to calculate the tSZ 
power spectrum are the halo mass function and the pro- 
jected radial electron pressure proflles for clusters over a 
wide range of mass and redshift. An important test is 
therefore to compare directly the three-dimensional pres- 
sure proflles of our model against those measured from 
X-ray observations. 

Arnaud et al. (2010, henceforth AlO) measured the 
electron pressure profile for intermediate to high-mass 
(10^'' < M§^1^/Mq < 10^5), low redshift {z < 0.2) clus- 
ters in the REXCESS sample (Bohringer et al. 2007; 
Pratt et al. 2009). The scaling of the amplitude of 
the profile with cluster mass at i?5oo was found to be 
Pe oc j\f0.69±o.i6 g^j^j^ |.]-^^g jj-^ agreement with self-similar 

expectations (P500 oc M^/^, where P500 is the charac- 
teristic pressure in the self-similar model as described in 
Appendix A of AlO). Due to the low-redshift nature of 
the REXCESS sample, AlO were not able to constrain 
significantly the redshift evolution of the sample, assum- 
ing that it follows the self-similar form, F500 oc h{z)^^^. 
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Fig. 4. — Comparison between the observed and model M|qq^ — 
fg relation for increasing values of the feedback parameter ef = 
10-'^, lO-*^, 5 X 10-6, and IQ-^ (solid lines, from top to bot- 
tom). Points with error bars represent individual cluster obser- 
vations from Sun et al. (2009) (blue) and Vikhlinin et al. (2009) 
(red). The dashed line represents the minimal feedback model 
(tf = EDM = 0.0; i.e., star-formation only). The dot-dashed line is 
the maximal feedback model with ef = 10—"" and edm = 0.1. The 
horizontal dotted line represents the universal baryon fraction. 

Having rescaled the profiles to remove any mass depen- 
dence there was found to be little dispersion (approxi- 
mately 20%) around the mean profile outside of the core 
region, r > 0.2i?5oo. Within 0.2i?5oo, AlO observed a 
much larger dispersion in the measured pressure profiles, 
with the shape of the profile related to the dynamical 
state of the cluster. AlO found there to be a good agree- 
ment between the observed profiles and those measured 
from the hydrodynamical simulations of Borgani et al. 
(2004); Nagai et al. (2007a) and Piffaretti & Valdarnini 
(2008). As the observations of AlO extend only out to 
i?5oo, the simulation data were used to extend the best-fit 
profile out to 4i?5oo. 

3.3.1. Impact of Energy Feedback 

In the left panel of Figure 5, we compare the electron 
pressure profiles for our gas model over a range of val- 
ues of et against the AlO profile. We plot the pressure 
profiles scaled as Pe(?')/-P5oo(^/^50o)"^ to allow a clear 
comparison of the profiles in the outer regions, which 
contribute significantly to the tSZ power spectrum at 
the angular scales of interest. The solid blue line rep- 
resents the P09 profile within i?5oo, the radius within 
which it was observed. The dashed blue line represents 
the region in which the profile was determined from sim- 
ulations rather than observations. The shaded region 
denotes the 20% dispersion observed by AlO around the 
mean profile. We plot Pe(r) for two different masses, 
= 3 X 10^'^ h-^MQ (black dashed lines) and 
M^oo^ = 3 X 10" H-^Mq (red solid), both at z = 0.1. 
We note that the AlO profile was measured for clusters 
of mass M^^Q^ > 7 x 10" H-'^Mq and so the lower of the 
plotted masses represents an extrapolation of the mass 
dependence of this profile. As in Figure 4, we plot model 
profiles for ef = lO"'^, 10-^,5 x 10"^ and 10"^. ao is 



fixed at the fiducial value (0.18). 

For Mf§^ = 3 X H-^Mq (black dashed lines), 
varying the feedback parameter has a very small effect 
on the pressure profile, especially within i?5oo- As noted 
above, for higher mass clusters the feedback energy is 
a small fraction of the total binding energy and there- 
fore does not strongly infiuence the gas density and tem- 
perature distribution. This is not the case for lower 
mass clusters. For = 3 x 10^'^ h~^MQ, increasing 

the feedback lowers the overall electron pressure within 
2i?5oo, with the effect strongest in the central regions 
of the cluster. As demonstrated in Figure 4, increasing 
the feedback parameter has the effect of inflating the gas 
mass distribution, reducing the gas density within i?5oo 
and increasing it at larger radii. Although increasing the 
feedback parameter produces a small rise in gas tempera- 
ture, this is outweighed by the decrease in the gas density, 
thus lowering the overall thermal pressure within i?5oo- 
We note that at larger radii (r > 1.2i?5oo), the electron 
pressure increases with ef. 

For high-mass clusters we find that all values of ef pro- 
duce an excellent match to the AlO pressure profile in the 
range 0.1 < r/R^oo < 2. For M500 = 3 x 10" H'^Mq, a 
significantly lower mass than any of the observed clusters 
in the AlO sample, the amplitude of the pressure profile 
is consistently below that of the AlO profile within i?5oo, 
but above at larger radii. 

3.3.2. Impact of Non-thermal Pressure Support 

In the right panel of Figure 5, we show the depen- 
dence of the gas model pressure profile on the non- 
thermal pressure support parameter ao. We vary ao in 
range < ao < 0.30 in steps of 0.06 (from top to bot- 
tom at i?5oo)- As in the upper panel, the two sets of 
curves represent clusters of mass = 3 x lO^'^ and 

3 X 10^^ h~^MQ, respectively. The feedback parameter 
is fixed at lO^'', which we henceforth take as our fiducial 
value. 

As the fraction of non-thermal pressure support is in- 
creased, the profiles become steeper, significantly reduc- 
ing the thermal pressure in the outer regions. This 
is entirely expected from Equation 18. For M|J,q^ = 
3 X 10^^ h~^MQ, setting ao = 0.18 produces a good fit to 
the AlO profile all the way out to 3i?5oo. In general, all 
the values of ao explored produce a pressure profile that, 
for r < R500, lies within the 20% dispersion observed 
around the AlO profile. For M500 = 3 x 10" /i"^Mq the 
profiles tend to lie below the AlO profile within i?5oo, but 
predict a higher pressure beyond this radius. 

3.3.3. Mass and Redshift Dependence 

At r — i?5oo, our fiducial model (ef — 10^^ and ao — 
0.18) produces a mass scaling of Pe(i?50o) oc M° '^5±° oi. 
At 0.5i?5oo, this steepens to M° ™±°-°^ This steepening 
of the mass scaling at smaller radii was also observed by 
AlO, who found a scaling of A/O.eoio.ie ^,^o.78±o.i6 
at r = i?5oo and 0.5i?5oo, respectively^. 

In Figure 6, we plot the redshift evolution of the pres- 
sure profile where the sohd, dashed, dot-dashed, and dot- 
ted lines correspond to z = 0, 0.5, 1 and 2, respectively. 

® We assume the same error on the mass scaling at O.S-Rsoo as 
was measured at iJsoo for the AlO profile. 
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Fig. 5.— (Left) Electron pressure profiles for clusters of mass M|J|^ = 3x10^'' H-'^Mq (black dashed lines) and Mf^^ = 3x10^^ H-'^Mq 
(red lines), using feedback parameters values of ef = 10~^, 10~^, 5 X 10~®, and 10~^ (from bottom to top at 2iJ5oo). The blue line represents 
the AlO pressure profile, with the shaded region denoting the observed 20% dispersion around this profile. (Right) Same as the upper 
panel but for increasing values of the non-thermal pressure support parameter < ao < 0.30 in steps of 0.06 (from top to bottom) 



The blue lines represent the AlO profiles at each red- 
shift, for which the normalization is assumed to scale 
self-similarly, i.e., P500 oc /i(z)^-^''. 

At r « 1.5i?5oo our model reproduces the self-similar 
scaling of the AlO profile, but scales differently with red- 
shift at smaller and larger radii. As redshift increases, 
the profile becomes more centrally concentrated, with 
the pressure increasing with respect to the AlO profile 
for r < 1.5i?5oo, but decreasing at larger radii. At 
r = 0.5, 1 and 2i?5oo we find that the profile scales as 

oc /l(z)2-71±0.04^ ;j(^)2.79±0.05 ^nd /l(z)2-57±0.06_ 

There are three components of our fiducial model that 
evolve with redshift at fixed halo mass: the stellar mass 
fraction M^,/M5oo (and thus the feedback energy), halo 
concentration - both of which decrease towards increas- 
ing redshift - and the non-thermal pressure parameter 
a, which initially increases as (1 -I- z)°-^ but asymptotes 
toward an upper limit of a{z) = 0.33 at high redshift 
(Section 2.2.5). While the increasing non-thermal pres- 
sure support produces a negative redshift evolution, this 
is compensated to some extent by the decreasing concen- 
tration. At fixed mass, as halo concentration decreases 
so does the central gravitational potential, reducing the 
pressure in the inner regions as well as increasing it at 
larger radii. We find that removing the redshift evolu- 
tion of the concentration parameter (Equation 6) results 
in a significant weakening of the redshift evolution rela- 
tive to self-similar, Pe(2i?5oo) oc /i(z)^-^**='=°"^. However, 
the evolving stellar mass fraction does not have a strong 
effect on the redshift scaling of the pressure profiles. 

4. SUNYAEV-ZEL'DOVICH ANGULAR POWER SPECTRUM 

The primary goal of this work is to investigate the vari- 
ations in the predicted tSZ power spectrum caused by 
uncertainties in ICM physics, in this case parameterized 
by the energy feedback and non-thermal pressure support 
parameters described in Sections 2.2.4 and 2.2.5. In the 
previous section we demonstrated the effect of varying 
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Fig. 6. — Pressure profiles for model clusters at different red- 
shifts, where z = 0, 0.5, 1, and 2 (solid, dashed, dot-dashed, dot- 
ted and thin solid lines, respectively). The blue line is the ob- 
served pressure profile of AlO. The cluster mass is fixed at 



M5HSE = 3 X 1014 



IQ. 



these parameters on the pressure profile for individual 
clusters, comparing with results from high-quality X-ray 
cluster observations. We now repeat this exercise for the 
tSZ power spectrum. We also compare our fiducial model 
with previous simulations or models and the recent SPT 
observations of the tSZ power spectrum. 

4.1. Impact of Cluster Physics 

Figure 7 illustrates the tSZ power spectrum obtained 
while varying the energy feedback parameters (ef and 
cdm), the normalization of the halo mass-concentration 
relation (Ac), the non-thermal pressure support param- 
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Fig. 7. — (top, left) The tSZ power spectrum for different values of ef and cdm- At £ = 10,000, the values of ej are (from top to 
bottom) 10~'^,10~®,5 X 10~^ and 10~^. The dashed line denotes eoM = ef = 0. The dot-dashed line indicates a model with Cf = 10^^ 
and eoM = 0.1. (top, right) tSZ power spectrum for different normalizations of the halo concentration-mass relation. At £ = 10,000, 
Ac = 0.8, 0.9, 1.0, 1.1, 1.2. (from bottom to top) (bottom, left) tSZ power spectrum for increasing values of ag, where, from top to bottom, 
< ao < 0.30 in steps of 0.06. (bottom, right) Varying the value of the non-thermal pressure support redshift evolution parameter /3 in 
the range —1 < /3 < 1 (from top to bottom) in steps of 0.5 . In each case, the thick red line represents our fiducial model. 



Value 




EDM 


ao 


/3 


Ac 


Fiducial 


lO"'' 


0.05 


0.18 


0.5 


1.0 


Min. 











-1 


0.8 


Max. 


10^^ 


0.1 


0.3 


1 


1.2 



TABLE 1 
Range of Gas Model Parameters. 



eter (ap) and its redshift dependence {/3). In each plot, 
the thick red hne denotes the power spectrum obtained 
from our fiducial model parameters. Table 1 summarizes 
the fiducial, maximum and minimum of each parameter 
explored. The range within which we vary ao, Cf and edm 
are chosen such that our model brackets the observations 
described in the previous section. 

We plot the band powers in terms of = £{£ + 
l)Ci/i2n) in units of ^iK^ at 147 GHz (i.e., f{x^) = -1 
in Equation 2). Note that we truncate our model profiles 
at 2i?vir(~ 4i?5oo). This radius was chosen to roughly co- 
incide with that of the accretion shock that is observed 
in hydrodynamical simulations (e.g., Molnar et al. 2009). 
We do not expect our model to be valid beyond the shock 
radius. Furthermore, the contribution of the low-density 



gas outside of groups and clusters to the tSZ power spec- 
trum is expected to be small (Hallman et al. 2009; Trac 
et al. 2010). Nevertheless, by doubling (tripling) the 
outermost radius we find an average increase of only 
6% (7%) in tSZ power. 

In the top-left panel, we plot the tSZ power spec- 
trum varying the feedback parameter between ef = 
10"'^, 10"^, 5 X 10"^ and 10"^ (solid lines from top to 
bottom dX t — 10, 000). The dashed line represents ef — 
cdm = and the dot-dashed hne ef = 10~^,eDM = 0.1. 
Increasing ef and eDM has the effect of suppressing power 
at small angular scales, causing Vi to peak at increas- 
ingly large scales. For example, for ef = eDM = 0, Vi 
peaks at £ — 3000, while for the largest feedback val- 
ues probed, peaks a,t £ = 2000. Similar results have 
recently been reported by Battaglia et al. (2010), who 
investigate the impact of AGN feedback in hydrodynam- 
ical simulations on the SZ power spectrum. The suppres- 
sion of small-scale power is caused by the inflating effect 
of energy feedback flattening the pressure profiles in the 
inner regions of clusters. As demonstrated in Figure 5, 
increasing feedback has the strongest effect on lower mass 
clusters M < 10^"' H^^AIq, which principally contribute 
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to the power spectrum on small angular scales. 

In the lower-left panel of Figure 7, we show the tSZ 
power spectrum for increasing values of the non-thermal 
pressure parameter, ao- We increase ccq in steps of 0.06 
between and 0.30 (from top to bottom at ^ = 1000), 
while holding the other gas model parameters at their 
fiducial value. 

Increasing ao results in a significant decrease in power, 
particularly at large angular scales; at £ = 3000, the full 
range in a.o encompasses more than a factor of 2 in power. 
In Section 3.3 we demonstrated that increasing the non- 
thermal pressure parameter significantly decreases the 
thermal pressure at large radii. This drives the large 
reduction in tSZ power, particularly at low £. We note 
that for higher values of ao, a{z) rapidly reaches the 
maximum value of 0.33 for the fiducial value of /3 = 0.5. 
The effect of increasing ao therefore begins to saturate 
for ao > 0.18. This is particularly noticeable on small 
scales {£ > 5000), where high redshift groups and clusters 
contribute significantly. It is evident that determining 
the level of non-thermal pressure support in groups and 
clusters, from all sources, is important in predicting the 
amplitude of their tSZ signal. 

In the lower-right panel of Figure 7, we plot (from 
top to bottom) the tSZ power spectrum for ^ = 
-1, -0.5, 0, 0.5 and 1. We hold ao fixed at 0.18. It is clear 
that the largest absolute change in power is at £ = 2500, 
although the fractional change is similar across the full 
range of angular scales plotted. This is in contrast to 
varying ao, which produces an increasing effect toward 
lower £. The reason for this is that varying /3 predomi- 
nantly affects high redshift objects, which contribute to 
the tSZ power spectrum at smaller angular scales than 
lower redshift objects of the same mass. However, ao pri- 
marily influences the pressure profiles at large radii (i.e., 
larger angiilar scales). Hence, there is a trade-off with 
the largest change in power being found at intermediate 
scales. 

In the upper-right panel we plot the tSZ power spec- 
trum while varying the normalization of the mass- 
concentration relation around the Duffy et al. (2008) 
value. This is achieved by setting Ac = 0.8,0.9, 1.0, 1.1 
and 1.2 (from bottom to top at £ = 10, 000), where Ac is 
defined in Equation 6. The results show that increasing 
halo concentration boosts the SZ power at small angular 
scales. For a given mass, higher concentration halos have 
a deeper central potential, thus increasing the pressure 
required for the inner regions of clusters to maintain gas 
in hydrostatic equilibrium. This steepens the pressure 
profiles in the inner regions of clusters, resulting in the 
observed increase in power at small angular scales. It 
is interesting to note that the impact of varying Ac on 
the tSZ power spectrum is somewhat similar to that of 
varying the feedback parameter ef. 

4.2. Comparisons with Simulations 

In Figure 8, we compare our model to other recent 
models or simulations. When comparing with simula- 
tions, we plot both our fiducial model (dashed lines) and 
a second realization in which we adjust the model pa- 
rameters to match the physics assumed in each simula- 
tion (dot-dashed lines), for example, by turning off star 
formation and energy feedback. This enables us to exam- 
ine the impact of second-order cluster properties (such as 



morphology and substructure) that are not accounted for 
in our halo-model approach. 

In the top-left panel we compare our model with the 
simulations of Sehgal et al. (2010), who calculated the 
SZ power spectrum using simulated SZ sky maps gener- 
ated by applying the semi-analytic model of Bode et al. 

(2009) to halos identified in an N-body lightcone simula- 
tion. This approach thus includes the effects of variations 
in halo density profiles, scatter in the mass-concentration 
relation and halo morphology (although the requirement 
of hydrostatic equilibrium will suppress the impact of 
substructure). The key difference between the model of 
Bode et al. (2009) and that described here is our inclu- 
sion of radially dependent non-thermal pressure support. 
The dot-dashed line demonstrates that we are able to 
reproduce almost exactly the Sehgal et al. (2010) tSZ 
power spectrum when we set ao = 0. This implies that 
variations in dark matter halo structural properties, do 
not strongly affect the tSZ power spectrum at the angu- 
lar scales investigated here. Comparing the Sehgal et al. 

(2010) prediction with our fiducial model highlights the 
importance of incorporating non-thermal pressure sup- 
port in ICM gas models. 

In the top-right panel, we compare our fiducial power 
spectrum with that predicted by the model of Komatsu 
& Seljak (2002). This model also assumes that the ICM 
is in hydrostatic equilibrium, but does not include star 
formation or any non-gravitational heating mechanisms. 
While it is difficult to isolate the exact cause of the dif- 
ference between the two profiles, Battaglia et al. (2010) 
have recently demonstrated the Komatsu & Seljak (2002) 
model significantly over-predicts the thermal gas pres- 
sure at large radii in comparison to observations, which 
is consistent with the factor of 2 increase in power com- 
pared to our model. 

In the lower-left panel we compare our model with 
SZ power spectra measured from recent hydrodynamical 
simulations. The red-solid line shows the non-radiative, 
GADGET-2 simulation of Battagfia et al. (2010), and 
the red-dashed line shows the result of a rerun of this 
simiilation including cooling, star formation and AGN 
feedback. The dotted lines show the expected sample 
variance around the non-radiative simulation given the 
simulated map size (1.6° x 1.6°, see also Shaw et al. 
2009). The black dashed line shows our fiducial model. 
For comparison, we also plot our model prediction having 
set the stellar mass and feedback energy to zero (black, 
dash-dotted line). 

The simulation without star formation, cooling and 
feedback predicts a higher tSZ signal than the one that 
includes these processes. On small scales, this is princi- 
pally due to the suppression of the SZ signal in the cen- 
tral regions of lower mass structures as gas is expelled by 
AGNs (Battaglia et al. 2010). On larger scales the differ- 
ence is driven by the lower gas mass in clusters due to star 
formation. Our model without star formation or feed- 
back demonstrates reasonable agreement with the non- 
radiative simulation at intermediate scales {£ w 4000). 
However, at larger scales it lies systematically above the 
simulation, and below at smaller scales. The latter may 
be due to additional power in the simulation power spec- 
trum due to the presence of substructures. The difference 
at large scales is most likely due to an artificial suppres- 
sion in the abundance of the most massive clusters due 
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Fig. 8. — Comparison with simulations and other models (solid lines). In each panel, the black dashed line is our fiducial model and 
the dash-dot line our model having modified the input parameters to reflect the level of physics incorporated in the simulation (where 
applicable). In the lower- left panel, the red lines show the results of the hydrodynamical simulations of Battaglia et al. (2010). All results 
are plotted at 147 GHz and are scaled to erg = 0.8 and Q.m = 0.264. 



to the Umitcd simulation volume. 

We find that our fiducial model matches that of the 
AGN simulation at £ = 1800, but, as with the non- 
radiative simulation, produces more power at larger 
scales and less at smaller scales. We note that by increas- 
ing the amplitude of the mass-concentration relation by 
20% our model very nearly reproduces the Battaglia et al. 
(2010) simulations for i > 3000. The overall reduction 
between the cases with and without baryon cooling and 
feedback are similar for the model and simulations. This 
indicates that the feedback prescriptions incorporated in 
our fully analytic model are able to reproduce the results 
of sub-grid models of AGN feedback in (computationally 
intensive) hydrodynamical simulations. 

Finally, in the lower-right panel we compare our fidu- 
cial model with that predicted by the Arnaud et al. 
(2010) pressure profile. Our model lies below that in- 
ferred by the Arnaud et al. (2010) profile, particularly 
at small scales. As shown in Figure 5, group-mass ob- 
jects have a significantly lower pressure at small radii 
(r < -R500) in our model compared to the AlO profile, 
which results in less power at small angular scales. By 
reducing the level of non-thermal pressure support in our 
model by one third (i.e, = 0.12) and setting the red- 
shift evolution of this parameter to zero our model will 
produce a power spectrum very similar to that of the 



Arnaud et al. (2010) profile. 

4.3. Comparisons with Observations 

The current generation of SZ surveys (SPT, ACT and 
Planck) have the necessary combination of sensitivity, 
angular resolution, mapping speed and frequency cov- 
erage to make precise measurements of the tSZ power 
spectrum. Multifrcquency observations are important for 
measuring the tSZ amplitude as both radio point sources 
and dusty star-forming galaxies (DSFGs) provide signifi- 
cant foreground signals near 150 GHz. While the Poisson 
component of these point source populations has a dif- 
ferent angular dependence to the SZ signal and can thus 
be separated and removed, the distribution of DSFGs on 
the sky are expected to be clustered with a power spec- 
trum of similar shape to the SZ signal (Viero et al. 2009; 
Hall et al. 2010). Hence, for single frequency surveys, the 
clustered DSFG and SZ power spectra cannot be distin- 
guished. Fortunately, the spectral dependence of these 
signals differs greatly and thus observing at two frequen- 
cies allows the clustered DSFG and SZ components to 
be separated. With three frequencies the kinetic and 
tSZ signals can also be independently measured. 

In Figure 9, we compare our fiducial model of the tSZ 
power spectrum to the recent SPT measurements of the 
small-scale CMB power spectrum. In this figure, the blue 



Impact of Cluster Physics on the SZ Power Spectrum 



13 



points with error bars represent the DSFG-subtracted 
power spectrum measured by Lueker et al. (2010). This 
was obtained by combining the maps at 220 and 150 GHz 
to optimahy remove the point source signal. The result- 
ing DSFG-subtracted power spectrum has been rescaled 
(by a factor of 2.2) to preserve the amplitude of the pri- 
mary CMB signal (see Section 6.2 of Lueker et al. 2010). 
Therefore, in order to enable a comparison with the ob- 
servations we have multiplied the tSZ power spectrum 
calculated with our model by this same factor. The 
dashed black line is the best-fit CMB-HtSZ+kinetic SZ 
(kSZ) model (plus a residual point source component) 
to the SPT observations assuming no theoretical uncer- 
tainty in the SZ predictions. Lueker et al. (2010) adopted 
the simulations of Sehgal et al. (2010) to infer a value of 
CTg = 0.746 from this fit (i.e., from the amplitude of the 
tSZ power spectrum alone). The solid and dashed red 
lines are the tSZ power spectrum predictions for our fidu- 
cial gas model parameters for cts — 0.8 and ag = 0.775, 
respectively (keeping the primary CMB, kSZ and point 
source components fixed). As in Lueker et al. (2010), 
we use the 'homogeneous' kinetic SZ power spectrum of 
Sehgal et al. (2010). 

Figure 9 demonstrates that our fiducial model matches 
the best-fit SPT power spectrum (dashed black line) 
when we lower as to 0.775, reducing the discrapency with 
the WMAP inferred value and constraints from cluster 
abundances. Increasing the values of ao, ef or cdm would 
increase the value of erg derived from the SPT measure- 
ments further. In general, we find that our fiducial model 
scales as Ce oc erf '* (at £ = 3000), which is steeper than 
the power of seven scaling found by Komatsu & Seljak 
(2002). This is due to the larger fractional contribution 
of higher mass objects to the tSZ power at these scales in 
our model, the abundance of which is extremely sensitive 
to the value of erg- 

5. DISCUSSION AND CONCLUSION 

The primary goal of this paper is to investigate the 
level of theoretical uncertainty on the amplitude of the 
tSZ angular power spectrum due to uncertainties in the 
physics of the ICM. To this end, we have developed an 
analytic model for the pressure and density distribution 
of the ICM. Our model is based on that of Ostriker 
et al. (2005) which assumes the ICM resides in hydro- 
static equilibrium in the potential well of NFW halos 
with a polytropic equation of state. The model accounts 
for star formation, energy feedback and energy trans- 
fer from dark to gas during mergers. We have made a 
significant improvement to the model by accounting for 
non-thermal pressure support in the ICM due to random 
gas motions, calibrating the non-thermal pressure pro- 
file from hydrodynamical simulations. This method al- 
lows us to rapidly generate theoretical thermal SZ power 
spectra, enabling a thorough investigation of the relative 
impact of astrophysical and cosmological parameters and 
their degeneracies. 

Our model has four key free parameters: an energy 
feedback parameter ef, which determines the amount of 
non-gravitational heating of the ICM from supernovae 
or AGNs, a dark-matter energy transfer parameter eDM, 
which governs the energy exchanged from dark matter 
to gas during mergers, a non-thermal pressure support 
parameter, ao, which gives the level of non-thermal pres- 
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Fig. 9. — Comparisons with SPT measurements of the CMB 
power spectrum at small angular scales. The blue points with er- 
ror bars are the DSFG-removed power spectrum measured by SPT, 
multiplied by 2.2 to preserve the amplitude of the primary CMB 
power spectrum (see text). The dashed black line is the best-fit 
CMB-(-tSZ-(-kSZ model to the SPT results. The red solid line is 
our fiducial tSZ model (plus primary CMB, kSZ and residual point 
sources) with erg = 0.8, the red dashed line is our fiducial model 
with (Tg = 0.775 (almost directly on top of the black-dashed line). 
Results are at 153 GHz. The thin black line is WMAP5 best-fit 
primary CMB temperature anisotropy power spectrum. 

sure at z = and its redshift evolution. These param- 
eters were calibrated against low redshift X-ray obser- 
vations of groups and clusters, including the M500 — fg 
(gas fraction) relation (Vikhlinin et al. 2006; Sun et al. 
2009) and observed electron pressure profiles (Arnaud 
et al. 2010). These studies derive mass estimates for 
observed clusters using the equation of hydrostatic equi- 
librium, which simulations suggest underestimate true 
cluster mass by 10% — 20%. Therefore, to ensure a con- 
sistent comparison, we are careful to use the hydrostatic 
(rather than true) mass predicted by our model when 
comparing with observations. 

We have found that increasing the energy feedback pa- 
rameter, Cf, from 10~^ to 10 has little effect on the 
pressure and density of high-mass clusters (where feed- 
back energy is a small fraction of the total binding en- 
ergy), but produces a factor of 2 decrease in the gas frac- 
tion for group mass objects. Similarly, increasing the 
feedback energy dramatically reduces the gas pressure 
within i?5oo and increases the pressure at larger radii in 
groups. Raising the level of non-thermal pressure (i.e., 
increasing ao) steepens the pressure profiles significantly 
in the outer regions (r > i?5oo) of both groups and clus- 
ters. We found that a fiducial model of ef = 10~^ and 
ao — 0.18 provides an excellent fit to the thermal pres- 
sure profile observed by AlO. 

By combining our gas model with the halo mass func- 
tion, we investigated the impact of varying the energy 
feedback, non-thermal pressure support and its redshift 
evolution on the tSZ power spectrum. We also evaluated 
the effect of perturbing the normalization of the mass- 
concentration relation in our model around the mean re- 
lation reported by Duffy et al. (2008). We found that 
increasing ef (and eDM) suppresses power on small angu- 
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lar scales, but does not strongly affect large-scale power. 
Hence, the peak of the power spectrum shifts to larger 
angular scales as we increase the amount of energy feed- 
back. This is because ef has a more significant effect on 
lower mass systems, which primarily contribute power 
on small scales. Varying the normalization of the mass- 
concentration relation Ac has a somewhat similar effect; 
the power on small scales increases with Ac- Raising 
the concentration of the host halo deepens the central 
potential and thus steepens the gas pressure profile so as 
to maintain hydrostatic cqiiilibrium. 

Varying the non-thermal pressure parameter qq has 
the opposite effect to the feedback parameters; increas- 
ing Qo strongly suppresses power on intermediate and 
large scales {£ < 5000). The peak of the power spectrum 
shifts slightly to smaller scales as ao is increased. Gener- 
ally, we found that non-thermal pressure support, at the 
level observed in recent hydrodynamical simulations, has 
a significant effect on the tSZ power spectrum, reducing 
the amplitude by a factor of 2 relative to the thermal 
pressure-only case. 

We have demonstrated that our model reproduces the 
tSZ power spectrum measured from simulated maps con- 
structed by applying the semi-analytic model of Bode 
et al. (2009) to halos identified in an N-body lightcone 
simulation (Sehgal et al. 2010). This demonstrates that 
deviations in halo structural properties (such as scatter 
in the mass-concentration relation) do not strongly af- 
fect the tSZ power spectrum, justifying our halo model- 
based approach. However, comparing with hydrodynam- 
ical simulations demonstrates that analytic models may 
potentially miss some small scale power due to substruc- 
tures in simulations. A more detailed comparison of hy- 
drodynamical simulations with semi-analytic models is 
required to isolate the impact of second order cluster 
properties on the SZ power spectrum. 

Recent SPT measurements of the small-scale CMB 
temperature anisotropy power spectrum (Lueker et al. 
2010) have demonstrated there to be some tension be- 
tween the value of as implied by the measured ampli- 
tude of the tSZ power spectrum at ^ = 3000 and those 
derived from WMAP observations of the primary CMB 
power spectrum (Dunkley et al. 2009; Komatsu et al. 
2010), or from cluster abundances (Vikhlinin et al. 2009; 
Mantz et al. 2010). This tension can be alleviated if the 
simulations used by Lueker et al. (2010) overestimate the 
amplitude of the SZ power spectrum by approximately 
a factor of 2. When compared with the SPT results, we 
found that our fiducial model infers as = 0.775, and thus 
reduces the discrepancy in the value of as between this 
and other probes. Generally, we find that our fiducial 
model scales as Ce oc erf'*. 

The results presented by Lueker et al. (2010) encom- 
passes roughly 5% of the expected SPT final survey 
area. Over the next few years, SPT, ACT and Planck 
should produce precise measurements of the SZ power 
spectrum amplitude over a wide range of angular scales 
2000 <i<lO, 000. We have demonstrated that physics 
in cluster environments can modify the shape of the 
power spectrum, as well as the amplitude. By comparing 
the ratio of power measured at small, intermediate and 



larger scales it should be possible to derive information 
regarding the relative infiuence of feedback processes and 
non-thermal pressure support in the ICM (albeit over a 
wide range of mass and redshift). Given that a signifi- 
cant fraction of SZ power derives from high-redshift and 
low-mass objects, the tSZ power spectrum can provide 
an exciting tool for studying the state of gas in these 
objects. 

Prom a theoretical perspective, several improvements 

can be made to our modeling of the ICM. First, our pre- 
scription assumes feedback from supernovae and AGNs 
is even distributed into the ICM. Although a concerted 
effort to implement and study AGN feedback in hydro- 
dynamical simulations of cosmological volumes has only 
recently begun (Sijacki et al. 2007; Di Matteo et al. 2008; 
Booth & Schaye 2009; Teyssier et al. 2010), these studies 
will shed light on the extent to which non-gravitational 
heating sources are able to influence the thermal struc- 
ture of the ICM. While our model accounts for energy 
transfer between dark matter and gas during mergers, 
the efficiency with which the process occurs has been 
poorly studied in simulations. Furthermore, we have 
not accounted for any uncertainty in the slope and am- 
plitude of the stellar mass total mass relation used 
in our model, which is poorly constrained at higher 
redshifts. Finally, the evolution of the radial, mass 
and redshift dependence of non-thermal pressure support 
(from any source) needs to be studied in hydrodynam- 
ical simulations in more detail. Future X-ray observa- 
tories equipped with high-resolution calorimeters (such 
as ASTRO-H) will provide important constraints on the 
non-thermal pressure support due to internal gas motions 
via broadening of heavy ion emission lines. 

In this paper, we have not studied the degeneracies in 
our model between astrophysical and cosmological pa- 
rameters in the shape and amplitude of the tSZ power 
spectrum, and how these may limit the precision to which 
as, for example, can currently be constrained in this way. 
We leave this to a follow-up paper (Bhattacharya et al. 
2010b). However, our prescription allows the observed 
shape of the power spectrum to be used to weaken or 
break some of these degeneracies. Alternatively, one can 
marginalize over the feedback or non-thermal pressure 
parameters in our model in order to account for astro- 
physical uncertainty in the predicted tSZ power spec- 
trum. We intend to make the code for our model pub- 
licly available for use in analyzing observations of the 
small-scale CMB power spectrum. 
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